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Abstract. The irreversible growth of magnetic films is studied in three-dimensional confined geometries 
of size L x L x M, where M S> L is the growing direction. Competing surface magnetic fields, applied 
to opposite corners of the growing system, lead to the observation of a localization-delocalization (weakly 
rounded) transition of the interface between domains of up and down spins on the planes transverse to the 
growing direction. This effective transition is the precursor of a true far-from-equilibrium corner wetting 
transition that takes place in the thermodynamic limit. The phenomenon is characterized quantitatively by 
drawing a magnetic field-temperature phase diagram, firstly for a confined sample of finite size, and then 
by extrapolating results, obtained with samples of different size, to the thermodynamic limit. The results 
of this work are a nonequilibrium realization of analogous phenomena recently investigated in equilibrium 
systems, such as corner wetting transitions in the Ising model. 

PACS. 68. 08. Be Wetting - 68.35. Rh Phase transitions and critical phenomena - 05.70. Ln Nonequilibrium 
and irreversible thermodynamics - 75.70.Cn Magnetic properties of interfaces 

1 Introduction PP), while it is also an attractive, challenging phenomenon 

from the theoretical point of view. Indeed, many theoret- 
Certainly, the phenomenon of wetting is a problem of pri- i ca ] efforts involving different approaches have been de- 
mary importance in many fields of practical technological vo t e d to the study of wetting transitions at interfaces un- 
applications (such as lubrication, efficiency of detergents, der equilibrium conditions |2I3I4I5I6I7I8I9I10I11| (for gen- 
oil recovery in porous materials, stability of paint coat- era j overviews of the subject see e.g. |12I13] L and more 
ings, interaction of macromolecules with interfaces, etc. 
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recently some investigations on nonequilibrium wetting 2 Brief summary of recent studies on corner 

phenomena | 14 | 15 | 16 | have been carried out as well. wet ting under equilibrium. 

A related phenomenon, which has attracted great in- 

Since the present work addresses the phenomena of non- 

tercst in the past decade, is the so-called filling or cor- 

equilibrium corner wetting in a magnetic material, it is 

ner wetting transition of a fluid adsorbed on a wedge 

useful to briefly describe recent progress in the study of 

17 18 19 . At present, it is experimentally feasible to con- 

their equilibrium counterpart )2ll2^l2^l2^l2^IM27l2^9l3()) . 

trol the shape of solid surfaces at a nanoscopic level, and 

Figure ^ shows a sketch of the corner geometry, of size 

indeed a growing variety of new types of interfacial phase 

L x L, used for the study of corner wetting in two di- 

transitions have been found in fluids confined by such 



.yflns.L Of course, the same sketch can be considered 

and 3D |27I28| wedges has also recently been investigated 



structured substrates gD|. The critical wetting in 2D |21l22l2j0^6i 

as a transverse view of a three dimensional array of size 



L X L X M (with M 3> L), as used in the present work. 

theoretically under equilibrium conditions. In particular, 

Let us consider a wedge in contact with the gas phase 

it should be mentioned that the exact solution for a two- 

of a fluid at temperature T and chemical potential ji. 

dimensional rectangular fsing ferromagnet forming a cor- 

Based on thermodynamic (macroscopic) arguments 1221 . 

ner with a surface field applied to the spins on edges has re- 

__ . the theory of corner wetting predicts that the filling of 

cently been obtained |26j . Furthermore, many theoretical 

the wedge by the fluid occurs at the transition tempera- 
results have later been corroborated by extensive Monte 

_ _ ture TV, given by 122123) 

Carlo simulations performed in both 2D 29 and 3D 130) . 

Within the context of these recent developments, the 

7r (T / )= 7 r/2-0, (f) 

aim of this work is to investigate the phenomenon of corner 

wetting under far-from-equilibrium conditions by means where n p2] is the contact angle describing the droplets 

of a Monte Carlo approach. To our best knowledge, this on a planar interface in the regime of incomplete wetting 

paper presents the first study of nonequilibrium wedge and cj) is the angle between the wall and the diagonal of the 

wetting. It is also worth mentioning that the irreversible wedge (see Figure QJ. Consequently, T/ is smaller than the 

filling of cavities with magnetic materials is also a frontier wetting temperature T w characterising the liquid-vapour 

topic in the field of the physics of new materials, since interaction on a planar substratum. 

modern trends in technology |31j require the character- Similar arguments can be drawn for a magnetic sys- 

ization of the filling of templates, containing imprinted tern in a wedge geometry, such as the Ising magnet in the 

nanometer/micrometer sized features, with a depositing square lattice in the presence of competing confinement 

material [32133) . fields as shown in Figure ^ For the sake of simplicity, 
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let us discuss the filling transition (also known as corner- 
wetting transition) in d = 2 considering the case <j) — tt/A 
only. The transition is of second-order and occurs at a 
critical field H cw |26| . which depends on T, and it is ob- 
served at temperatures below the order-disorder Onsager 
critical temperature of the Ising magnet. At low temper- 
atures, the interface remains bounded to one wall and its 
mean distance to the corner ((l ), see Figure ^) diverges 
according to a power-law behavior when approaching the 
transition. Of course, the growth of the correlation lengths 
is bound by the finite lattice size and true divergences are Fig- x _ sketch of the simulation geometry where a transverse 




only possible in the thermodynamic limit. So, one has 



(l a ) oc t 



t = H CW (T) — H, 



(2) 



where /3 S = 1 is the order parameter critical exponent. The 
interface also develops correlations along the directions 
parallel and perpendicular to it (see Figure QJ. The cor- 
responding correlation lengths (t; x and respectively), 
also diverge at the critical point according to 



plane (k = constant) of size L x L of the square lattice is 
shown. Free boundary conditions are used for the spins in rows 
i — l,i — L and j = l,j = L, with integers £ [1,L] la- 

belling the lattice sites. Additionally, a magnetic field H acts 
on all spins in the rows i = 1 and j = 1, while a magnetic field 
— H acts on all spins in the rows i = L and j = L. A typical 
configuration of the system below the critical point then con- 
tains an interface between a domain of positive magnetization 
(+) and a domain with negative magnetization (— ). Fluctua- 
tions of the interface are characterised by correlation lengths 
(3) 6e an d £± in the directions parallel and perpendicular to the 
interface, respectively. 



where the correlation length exponents are v x = 1 and 
v±_ = 1, respectively. 

On the other hand, taking first the limit L — > oo and 
then the limit t —> + , the probability distribution of the 
mean position of the interface is expected to decay expo- 
nentially, namely 



(4) 



For finite systems of side L, equation J3J needs to be sym- 
metrized with respect to both corners to which the in- 
terface can be bound |2H]. Equation holds within the 
regime of incomplete corner wetting (or equivalently to 
the case of incomplete corner filling, when the interface of 
the magnetic domain filling the sample is still bound to 
one corner far from the diagonal crossing the sample from 
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left-top to right-bottom Q ). However, in the regime of complete wetting the interface is no longer bound to the 
complete corner wetting (or equivalently complete filling), wall (or corner, respectively). 
t < 0, the average location of the interface lies along the 
diagonal taken from the upper-left corner to the lower 
right corner. So, one has (l a ) = Lj\[2 and the probability 
distribution becomes 



3 The Magnetic Eden Model in a corner 
geometry and the simulation method 



For definiteness, we will adopt a magnetic language through- 

P((l )) ex / ro ~ (jo)] \ fr\ 

" °'' ^ ^ 2^ out, although the relevant physical concepts discussed here 

with [2"H| can ^ e ex t en ded to other systems such as fluids, polymers, 

and binary mixtures. The irreversible growth of a ferro- 



a«£», (6) 



magnetic material has been studied by using the so-called 
magnetic Eden model (MEM) an extension of the 

since it can be interpreted in terms of a random walk de- classical Eden mode l 37 m w hich the growing particles 
scription of the interface in d = 2, such that excursions in have an additional degree of freedom, i.e., the spin. Sev- 
the direction perpendicular to the average interface orien- eral recent investigations based on this modcl showed a 
tation add up randomly [33]. ricn var j e ty f interesting phenomena, such as the occur- 

Equation © implies that for complete wetting, the rence of morphological phase transitions associated with 
average magnetization vanishes ((to) = 0) and the Gaus- thc gr0 wing interface (SHJ, Ising-like order-disorder phase 
sian distribution of the interface location translates into a transitions spontaneous magnetization reversals 
Gaussian distribution of thc magnetization given by and a f ar -from-cquilibrium wetting transition driven by 

competing surface magnetic fields [15110] . 

/ to 2 1? \ 

P(m) ~ exp I , (7) The MEM in (2 + 1)— dimensions is studied in thc 



2k B T X L 

square lattice by using a rectangular geometry Lx Lx M 

with 

(with M ^> L). The location of each spin on the lattice 

XL « e± « L, (8) 

is specified through its coordinates (i,j,k), (1 < i, j < L, 
where the size-dependent width of the Gaussian is given 1 < k < M). A transverse plane of the used geometry, 
by xl = ((to 2 ) — (to) 2 ), with (m) 2 = within the complete taken for k = constant, is sketched in Figure^ The start- 
wetting regime. It should be noticed that, as in the case ing seed for the growing cluster is a plane of L x L parallel- 
of the localization-delocalization transition observed for oriented spins placed at k — 1, and cluster growth takes 
the Ising magnet in a slit geometry |2I44| . in the case of place along the positive longitudinal direction (i.e., k > 2). 
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Open boundary conditions are adopted in both transverse Since the observables of interest (e.g. the mean trans- 
directions, and competing surface magnetic fields H > verse magnetization and its fluctuations, see below) re- 
and W = —H are applied to opposite corners, as shown in quire the growth of samples with a large number of trans- 
Figure ^ Then, clusters are grown by selectively adding verse planes of size L x L, clusters having up to 10 9 spins 
spins (S^k — ±1) to perimeter sites, which arc defined have typically been grown for lattice sizes in the range 
as the nearest- neighbour (NN) empty sites of the already 12 < L < 96. It should be mentioned that the grow- 
occupied ones. Considering a ferromagnetic interaction of ing front leaves voids behind that are incorporated to the 
strength J > between NN spins, the energy E of a given bulk of the sample during some transient period. How- 
configuration of spins is given by ever, since these voids are also perimeter sites they are 

ultimately filled in during the growth process. Hence, far 



E = -- I s wkS, 



behind the active growth interface, the system is compact 



°v k °i j'k I and frozen. So, the different quantities of interest are mea- 

(ijk.i' j' k') 

(\ / \ sured on defect-free transverse planes. 
£ Sljk + Y, Silk J + H S L] k + S iLk (9) 
(jk) (zk) J \(jk) (zk) J 

where the summation (ijk,i j k ) is taken over occupied Moreover, the update algorithm is quite time consum- 
NN sites, while the remaining terms are sums over occu- ing, as compared to standard Ising simulations, since the 
pied surface sites, in order to take into account the effect growing probability has to be computed after each depo- 
of the surface magnetic fields. The Boltzmann constant is sition event. Hence, it should be noticed that the results 
set equal to unity throughout, and the temperature, en- obtained in this work involve a large computational effort, 
ergy, and magnetic fields are measured in units of J. The Let us also remark that, although both the interaction en- 
probability for a perimeter site to be occupied by a spin is ergy and the Boltzmann probability distribution consid- 
taken as proportional to the Boltzmann factor exp(— -^P), ered for the MEM are similar to those used for the Ising 
where AE is the change of energy involved in the addi- model with surface magnetic fields, these two models op- 
tion of the spin. At each step, the probabilities of adding erate under extremely different conditions, since the MEM 
up and down spins to a given site have to be evaluated describes the irreversible growth of a magnetic material, 
for all perimeter sites. After proper normalization of the while the Ising model deals with a magnet under equilib- 
probabilities, the growing site and the orientation of the rium. Previous studies have demonstrated that the MEM 
spin are determined with Monte Carlo techniques. For ad- in (1 + 1)— dimensions is not critical but exhibits a second- 
ditional details on the MEM and the simulation method order transition at T c = 0.69 ±0.01 in (2 + 1)— dimensions 
see e.g. [IBIITSI3tildHl4()| . EH- 
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4 Results and discussion 

The growth of magnetic Eden aggregates in a (eJ+1)— confii 
geometry is characterized by an initial transient period 
followed by a nonequilibrium stationary state that is in- 
dependent of the initial configuration ^S]. The length of 
the transient is of the order of L d (d = 2 in the present 
study), and the proportionality constant depends on T. 
We have checked that disregarding the first 100L 2 trans- 
verse planes, subsequent averages are independent of the 
initial configuration. So, all results reported in this paper 
are obtained under these conditions. 

Depending on the temperature and the surface mag- 
netic field considered, different regimes are clearly distin- 
guished upon growth of a magnetic aggregate in the con- 
fined corner geometry, as illustrated by the typical snap- 
shot configurations of transverse [k = constant) planes 
shown in Figure [21 Furthermore, in order to characterize 
these different phases quantitatively, we will locate the 
boundaries between them in a magnetic field versus tem- 
perature (H — T) phase diagram (see Figure^. Below, we 
will first describe the results for a confined, finite-size sys- 
tem, and then we will obtain the phase diagram in the 
thermodynamic limit by means of extrapolation proce- 
dures. 

Let us now discuss the localization-delocalization be- 
havior of the interface observed when the value of the 
surface field is fixed (H = 0.4) but the temperature is 
increased. For this qualitative description we will use the 
typical snapshot configurations shown in Figure |21 Figure 
EJa) shows the regime characteristic of low temperatures 



(a) (b) 




(c) (d) 




Fig. 2. Typical snapshot configurations of transverse (k — 
constant) planes. Black points correspond to down spins, while 
up spins are left in white. The surface field is positive (negative) 
on the left and bottom (right and top) sides of the wedge, as 
shown in Figure The snapshots correspond to lattices of size 
L = 96 and are obtained keeping the surface fields constant 
and varying the temperature, as follows: (a) H = 0.4, T = 
0.6; (b) H = 0.4, T = 0.65, (c) H = 0.4, T = 0.7, and (d) 
H = 0.4, T = 1.0. 

(i.e. far below the critical temperature of the MEM). Here, 
the interface between up and down domains is localized 
close to the upper-right corner. Of course, during very 
long-time simulations one observes the localization of the 
interface close to both corners with the same probability. 
Keeping the same field but increasing the temperature, 
the interface starts to depart from the corners and, close 
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to the effective corner wetting transition temperature of 
the considered lattice, it remains still localized close to the 
diagonal of the sample, as shown in Figure^b). However, 
further increase of the temperature causes the interface to 
become delocalized (see Figure Etc)). Here the domain of 
positive spins makes frequent excursions up to the bound- 
ary with negative field (see the lower-right corner of Figure 
[21c)). Also, notice that a relatively weaker excursion of the 
negative domain to the positive boundary field is observed 
at the upper-left corner. Of course, due to the symmetry Fi S- 3 - The H — T phase diagram corresponding to a lat- 




o.o 

0.5 



0.5 0.6 0.7 0.8 0.9 
T 



| corner wetting (this work) 
wetting (Ref. [15]) 



of the system, these excursions cause the temporal aver- 
age position of the interface to be just at the diagonal of 
the sample with (^ ) = L/\/2, while its probability distri- 
bution is given by equation Q. For this regime the total 
magnetization vanishes and its probability distribution is 
a Gaussian centered around m = 0, see equation J7|). Fi- 
nally, at higher temperatures and within the disordered 
phase (i.e., above the critical temperature of the MEM), 
the interface between domains of up and down spins van- 
ishes due to the thermal noise, as shown by the typical 
high-temperature configuration of Figure |2Jd) . 

Let us now discuss the interface localization-delocalization 
"phase diagram" of a confined system of lateral size L — 
12, as shown in Figure |21 

In order to locate the transition curve we have mea- 
sured the magnetization fluctuations, given by 



tice of size L — 12. The vertical straight line at T C (L) = 0.84 
shows the L— dependent critical temperature of the MEM in 
the absence of any magnetic field, which separates the low- 
temperature ordered phase from the high-temperature dis- 
ordered phase. Besides, the corner localization-delocalization 
transition curve (full circles and solid line) obtained in this 
work, as well as the localization-delocalization transition curve 
(dashed line) that corresponds to the MEM with competing 
surface magnetic fields applied to parallel confinement walls 
|15j . is shown for comparison. The inset shows the phase dia- 
gram corresponding to the thermodynamic limit, as obtained 
by extrapolating results obtained for systems of different lat- 
tice size. Again, the corner wetting transition (solid line) is 
compared to the wetting transition of the MEM (dashed line) 
obtained in Reference |15| . More details in the text. 



by 



X (L,T) = -((m 2 )-(H) 2 ) , (10) 

where (...) means the average taken over a large number 
of transverse planes in the stationary regime, and m = 



m 



1 L 

72 s v k 



L 2 



(11) 



It should be noticed that under equilibrium condi- 
tions and provided that the fluctuation-dissipation the- 



m(k, L,T, H) is the mean transverse magnetization given orem holds, the fluctuations of the magnetization can be 
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Ql 1 1 1 , 1 , 1 

0.6 0.8 1 1.2 

T 

Fig. 4. Plots of the magnetization fluctuations as a function of 
the temperature obtained for different values of the magnetic 
field, as indicated. Following a standard procedure |44|. the 
maxima of x(L, T) define the points of the corner localization- 
delocalization transition curve on the H — T plane, as shown 
in Figure |3 

identified with the susceptibility. However, this may not 
be the case for the present far-from-equilibrium system. 

Working under equilibrium and far from it 42 43 
is usual to locate L— dependent 'critical' points at the 
maximum of x(-^> T). This procedure has proved to be also 
useful to locate L— dependent localization-delocalization 
transition points (namely wetting "critical" temperatures 
T CW (L) ) under equilibrium |44| and far from it |15j . So, 
we have also used the same method in the present case. In 
order to illustrate this procedure, Figure 0] shows plots 
of \{L,T) as a function of the temperature for differ- 
ent values of the surface magnetic field. The maximum of 
each curve defines a point in the transition curve shown 
in Figure It should be noticed that the localization- 
delocalization transition in a confined system actually cor- 



-from-equilibrium magnetic growth model 

responds to a so-called quasi-corner wetting transition, 
which is indeed the precursor of the true corner wetting 
transition that occurs in the thermodynamic limit. 

In Figure|3| the vertical straight line at T C (L = 12) = 
0.84 shows the location of the L— dependent "critical" 
temperature of the finite system in the absence of surface 
fields. This value was obtained previously for the MEM in 
the absence of an external magnetic field and corresponds 
to the slit geometry [25]. Although, as it is well known 
from finite-size scaling theory, there is some degree of ar- 
bitrariness in locating the L— dependent critical tempera- 
ture T C (L) of a finite system, the actual critical point T c of 
the infinite system, which can be obtained by extrapolat- 
ing T C (L) to the L — > oo limit, is unique and independent 
of any particular choice for the finite-size critical point. Let 
us also recall that for low temperatures (T < T C {L)) the 
phase diagram corresponds to the ordered growth regime, 
while high temperatures (T > T C {L)) are associated with 
the disordered growth regime. 

For the sake of comparison, the localization-delocalization 
transition curve obtained in this work for the corner geom- 
etry is compared to the localization-delocalization transi- 
tion curve of the MEM (dashed line in Figure |3J) that 
corresponds to competing surface magnetic fields applied 
to parallel confinement walls i- e -> the so-called slit ge- 
ometry. Notice that in this case the critical curve ends, for 
H = 0, at the L— dependent critical point T C (L = 12) = 
0.84. It is observed that the phase diagrams obtained here 
are similar to the analogous phenomena of wetting and 
corner wetting observed in 2D equilibrium systems, as e.g. 
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in the confined Ising model |29I44| . In fact, for a given sur- 
face field, in the case of corner wetting, the correspond- 
ing critical temperature is always smaller than for pure 
wetting except, of course, for H — 0. This scenario is in 
agreement with equation Q), which was obtained on the 
basis of thermodynamic (equilibrium) considerations |22) . 

It should be noticed that the transition curve corre- 
sponding to the corner geometry intercepts the horizontal 
axis (H = 0) close to T C (L = 12) ~ 0.81, i.e., a value 
slightly smaller than the L— dependent order-disorder crit- 
ical temperature of the MEM model in the absence of sur- 
face fields. This shift simply reflects the different boundary 
conditions that one has to use to study both systems. In 
fact, for the corner geometry one has to apply open bound- 
ary conditions, while for studying the confinement due to 
parallel walls one has to adopt open (periodic) boundary 
conditions along (in the direction perpendicular to) the 
walls. Of course, both localization-delocalization transi- 
tions lie within the ordered growth regime where an inter- 
face between domains is well defined. 

Further insights can be gained by examining the be- 
havior of the probability distribution of the mean trans- 
verse magnetization, Pl{tti). In fact, in the thermody- 
namic limit, the order parameter probability distribution 
of an equilibrium system at criticality is universal (up 
to rescaling of the order parameter), and hence it con- 
tains information on the universality class of the system 
[45146147] . Figure [5] shows the dependence of Pl(to) on 
the surface magnetic field for fixed values of the temper- 
ature (T = 0.65) and the lattice size (L = 12). For larger 




Fig. 5. Order parameter probability distributions obtained for 
T = 0.65 and different values of the field, as indicated. The 
lattice size is L = 12. 

fields (H — 1.8 in FigureEJ) the distribution is a Gaussian 
centered at m = 0, in agreement with the fact that the 
average location of the interface (for temperatures below 
T C (L)) is a straight line crossing the sample from one cor- 
ner, where the magnetic fields have opposite direction, to 
the other (see e.g. the snapshots shown in Figures |21 (b) 
and (c), as well as the corresponding discussion). Hence, 
the average magnetization (m) = 0, and the Gaussian dis- 
tribution of the interface location given by equation |5| 
translates into a Gaussian distribution of the magnetiza- 
tion given by equation J7J). So, within this wet regime the 
width of the Gaussian can easily be measured and the re- 
lationship xl oc L (see equations J7J) and © ) has been 
verified, as shown in Figure [5J This behavior is similar to 
that observed in the interface localization-delocalization 
transition under equilibrium conditions, in both the slit 
and corner geometries 29 . 
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Fig. 6. Plot of xl versus L, as obtained from the order parame- Fig. 7. Extrapolation of the results to the thermodynamic 
ter probability distribution measured within the wet phase, for limit: (a) plot of T CW {L) versus L~ , for H cw = 0.2, and (b) 
H — 0.8 and T = 0.65. plot of H CW (L) versus L~ , for T cw — 0.5. The extrapolated 

fits (which take into account the larger lattices only, i.e., 32 < 



Decreasing the magnitude of the field, the onset of two 



L < 96) are also shown. 



distinct maxima shows gradually up (e.g. for H — 1.8 

in Figure 0), as is also observed in equilibrium systems m the general context of continuous phase transitions, due 

|2|29|M|44|. Indeed, this behavior not only reflects the fact to finite-size effects the effective location of second-order 

that the interface is located close to one of both corners transitions becomes shifted according to gT] 
with the same probability, but also it is the signature of 
a thermal continuous phase transition taking place at fi- 

T c (L)=T c (oo)+AL- 1 ^, (12) 

nite critical temperature [45146147] . Finally, well inside the 

nonwet phase, for H = 0.5 in FigureEJ), one has an expo- where T C (L) is the effective L— dependent "critical" tem- 

nential decay of the distribution as already reported for perature, T c (oo) is the actual critical temperature in the 

the equilibrium counterpart thermodynamic limit, A is a constant and v is the correla- 

Let us now consider the extrapolation of the L-dependent tion length exponent related to the proper order paramc- 

critical points to the L — > oo limit, in order to show that ter. In the present case, the location of the interface ((/ ))j 

the above discussed behavior of the interface is not only see equation (j5J), is the appropriated order parameter for 

characteristic of confined systems, but it is also present in the corner wetting transition, rather than the magnetiza- 

the thermodynamic limit, hence leading to the phase dia- tion as in the standard Ising (or MEM) model. The asso- 

gram shown in the inset of Figure|21 It is well known that, ciated correlation length is given by £j_ (see equation Ipty) 
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describing the growth of correlations in the direction per- 
pendicular to the interface (see Figure^). For the confined 
Ising system in the corner geometry under equilibrium one 
has that the predicted exponent is z^j_ = 1. So, we have 
performed the extrapolations to the thermodynamic limit 
by taking v = 1. Of course, an accurate numerical de- 
termination of the exponent is beyond our computational 
capabilities. 

So, let us now apply these concepts to the corner wet- 
ting transition of the MEM. Considering the L— depen- 
dent critical points (H CW (L), T CW (L)) on the localization- 
delocalization transition curves corresponding to finite sys- 
tems, we fixed the value of either the magnetic field or the 
temperature, and extrapolated the results to the L^ 1 — > 
limit. For the purpose of illustration, Figure EJa) shows 
a plot of T CW (L) versus L" 1 obtained for a fixed value of 
the magnetic field (H cw = 0.2), while Figure[7Jb) shows a 
plot of H CW (L) versus L~ x for a fixed value of the tempera- 
ture (T cw — 0.5). Considering only larger lattices (namely, 
within the range 32 < L < 96), the extrapolations allow 
us to determine the corner wetting transition curve cor- 
responding to the thermodynamic (L ~ * oo) limit, which 
is shown in Figure [3] Deviation from the linear behav- 
ior observed for smaller lattices may be due to finite-size 
corrections or eventually to the fact that the correlation 
length exponent could slightly depart from unity, as al- 
ready discussed above. It should be noticed that, as in the 
case of confined (finite-size) systems, the corner wetting 
transition (lower curve, solid line) is here compared to the 
wetting transition of the MEM (upper curve, dashed line) 
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obtained in Reference JS]- As observed in the confined 
Ising model |29I44| . the critical curve in the case of corner 
wetting lies also below than that of pure wetting, except 
for H = 0. 

After the evaluation of the critical points it is nat- 
ural to discuss the trends of the magnetization fluctua- 
tions upon varying the system size and, subsequently, an- 
alyze whether a scaling plot could be obtained. Figurc|^a) 
shows plots of x(L, T) as a function of the temperature ob- 
tained for different values of L and keeping constant the 
surface magnetic field H = 0.2. As expected, the different 
curves exhibit rounding and shifting effects characteristic 
of finite systems. For the corner wetting transition under 
equilibrium conditions it is known that the susceptibility 
of the finite square scales according to |2H] 

x(L,T)<x^x*(Lt), (13) 

where \* is a scaling function and t = T — T cw (see also 
equation (|10fl ). So, in the inset of figure|HJa) we show that 
the maxima of the magnetization fluctuations Xmax (mea- 
sured at T CW (L) ) of the non-equilibrium counterpart also 
scales with L 2 . Furthermore, figure |S{b) shows a scaling 
plot of the magnetization fluctuations according to equa- 
tion (|13fl . The data shown correspond to T > T cw and 
the expected asymptotic behavior (dashed-dotted line) is 
roughly achieved for the larger samples only. Of course, 
a data collapse of better quality would be desirable but, 
regrettably, it is beyond our computational capabilities. 
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rial under far-from-equilibrium conditions, using exten- 
sive Monte Carlo simulations. The occurrence of this phe- 
nomenon was firstly described at a qualitative level, by 
means of typical snapshot configurations that depend on 
the control parameters (surface magnetic field and tem- 
perature) , and subsequently quantitatively, by identifying 
well defined regions in the H — T phase diagram. After 
presenting the interface localization-delocalization transi- 
tion phenomenon in (finite-size) confined geometries, the 
results were extrapolated to the thermodynamic limit. 
These results, obtained in the framework of nonequilib- 
rium growth systems, are a novel realization of analogous 
phenomena, which have recently been observed in equilib- 
rium systems. Hence, we hope that this work will con- 
tribute to the understanding of the irreversible growth 
of binary mixtures in confined geometries, as well as of 
wetting-related phenomena. 
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